Wave packet dynamics in a monolayer graphene 
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The dynamics of charge particles described by Gaussian wave packet in monolayer graphene is 
studied analytically and numerically. We demonstrate that the shape of wave packet at arbitrary 
time depends on correlation between the initial electron amplitudes 0) and 7/;2(r, 0) on the 

sublattices A and B correspondingly (i.e. pseudospin polarization). For the transverse pseudospin 
polarization the motion of the center of wave packet occurs in the direction perpendicular to the 
average momentum pb = hko. Moreover, in this case the initial wave packet splits into two parts 
moving with opposite velocities along po . If the initial direction of pseudospin coincides with average 
momentum the splitting is absent and the center of wave packet is displaced at t > along the same 
direction. The results of our calculations show that all types of motion experience zitterbewegung. 
Besides, depending on initial polarization the velocity of the packet center may have the constant 
component Vc = uf{a), where u ~ 10* cm/s is the Fermi velocity and f{a) is a function of the 
parameter a — kod {d is the initial width of wave packet). As a result, the direction of the packet 
motion is determined not only by the orientation of the average momentum, but mainly by the phase 
difference between the up- and low- components of the wave functions. Similar peculiarities of the 
dynamics of 2D electron wave packet connected with initial spin polarization should take place in 
the semiconductor quantum well under the influence of the Rashba spin-orbit coupling. 

PACS numbers: 73.22.-f, 73.63.Fg, 78.67.Ch, 03.65.Pm 
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I. INTRODUCTION 



At last years the dynamics of wave packets in 
2D electron gas and other systems in solids in- 
cluding the phenomenon of zitterbewegung (ZB) or 
trembling motion has been the subject of numerous 
studiesJ-i^i^iii^i^ii^i^iiSiiiii^i^ Firstly the oscillatory mo- 
tion analogous to the relativistic zitterbewegung in two- 
dimensional systems with the structural and bulk inver- 
sion asymmetry was investigated by Schliemann et al.^ 
In the recent work by authors the detailed studying of 
the electron wave packet dynamics in the semiconductor 
quantum well under the influence of the Rashba spin- 
orbit coupling was performedii^ It was shown (analyti- 
cally and numerically) that the initial wave packet splits 
into two parts with different spin polarizations propagat- 
ing with unequal group velocities. It was demonstrated 
also that the splitting and overlapping of wave packets 
leads to the damping of zitterbewegung. 

As well known, the electron zitterbewegung in relativis- 
tic physics at first time was predicted by Schrodinger^^ 
(see alsoi^). This phenomenon is caused by the inter- 
ference between positive and negative energy states in 
the wave packet. The frequency of ZB motion is deter- 
mined by the gap between these two states and when the 
amplitude of oscillations in a particle position is of the 
order of the Compton wave length. This phenomenon 
was discussed also in Refs. [17-19]. 

In the papers by Rusin and Zawadzkii2ii^ the evolution 
of the wave packet in a monolayer and bilayer graphene 
as well as in carbon nanotubes was analyzed. The exact 
analytical expressions for two components of wave func- 
tion and average value of position operator were found 



for bilayer graphene, which allowed to obtain analytical 
results for the ZB of Gaussian wave packet. It was shown 
that the transient character ZB in bilayer graphene is due 
to the fact that wave subpackets related to positive and 
negative electron energies move in opposite directions, so 
their overlap diminishes with time. At the same time the 
dynamics of the wave packets in a monolayer graphene 
in Ref.[12] was not investigated fully. 

In this work we present the detailed description of 
wave packet evolution including the phenomenon of ZB 
of the packet center in a monolayer graphene. We inves- 
tigate the influence of the initial pseudospin polarization 
on the dynamical characteristics of the moving electron 
wave packet. Our analytically results are illustrated by 
a graphic presentation. 

The outline of this paper is as follows. The expres- 
sion for the wave function at t > is found in Sec. II for 
the arbitrary initial pseudospin polarization. The time 
dynamics of Gaussian wave packet and, in particular, 
the ZB phenomenon are described in Sec. Ill for differ- 
ent correlations between the initial electron amplitudes 
on the sublattices A and B. We conclude with remarks 
in Sec.IV. 



II. BASIC EQUATIONS 

Graphene is a single layer of carbon atom densely 
packed in a honeycomb lattice. The two-dimensional 
Hamiltonian describing its band structure has the 
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where m « 10^ cm/s, p = (pxTPy) is the momentum 
operator defined with respect to the centre of the valley 
centered at the corner of the Brillouin zone with wave 
vector K. Pauli matrices ct^ operate in the space of the 
electron amplitude on two sites {A and B) in the unit cell 
of a hexagonal crystal. This internal degree of freedom 
plays a role of a pseudospin. The Dirac-like Hamilton H 
determines the linear dispersion relation 
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Here p = ^pl + Py, s = 1 for the electron in the con- 
duction band and s = — 1 for the valence band ("hole" 
branch of quasiparticles) . The corresponding eigenfunc- 
tions are given by 
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where coefficients ci and C2 determine the initial pseu- 
dospin polarization. We suppose that the packet width d 
is much greater than the lattice period and consequently 
il>{r^Q) is smooth enveloping function. We suppose also 
that the most of the states in valence band are unfilled, 
that corresponds to negative Fermi level located far from 
Dirac point (see also^). Substituting Eqs.(8a, 8b) in 
Eq.(4) and using the expressions (6) and (7) we obtain 



V'i(r,t) 



1 



^\C1V + \C2\ 



-.{ci(f>i{r,t) - C2(t>2{-x,y,t)), 
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with = E^dlPy-, 
p 

The time-evolution of an arbitrary initial state ipi'^, 0) 
in Shrodinger representation can be found with the help 
of Green's function Gfj_i,{f, r') 



= (c20i(r,t)-f ci02(r,i)), (10) 



where, for notational convenience, (t>i.2{^,t) denote the 
functions 
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where fi^v = 1,2 are matrix indices, corresponding to 
the upper and lower components of 'ip{f,t). These com- 
ponents are related to the probability of finding electron 
at the sites of the sublattices A and B correspondingly. 
The standard expression for Green's function is 



s=±l 

Using Eq.(3) for ipp^s-nirTt) we find 
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Using the cylindrical coordinates in Eqs.(ll), (12) and 
integrating over the angular variable, we have 
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Let us represent the initial wave function by Gaussian 
wave packet having the width d and nonvanishing average 
momentum poy = ^^o 



X y e 2 cos{qt)Jo{q\/r^ ~ — 2iay)qdq, (13) 
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FIG. 1: (Color on line). The electron probability density 
p{r,t) — \ipif + |^/'2p for initial wave packet determined by 
Eqs.(8a), (8b) with ci = 1 and C2 ~ for a = kod = 1.2 at 
the time t = 7 (in the units of d/u). 

where Jo{z), Ji{z) are Bessel functions. For the sake 
of convenience we introduce in Eqs.(13), (14) and every- 
where below the dimensionless variables, measuring the 
distance in the units of initial width of wave packet d and 
time in d/u units. Besides, instead of the wave vector ko 
we consider the parameter a = kod. 



III. ZITTERBEWEGUNG OF GAUSSIAN WAVE 
PACKET WITH DIFFERENT PSEUDOSPIN 
POLARIZATION 

Now we describe the time dynamics of Gaussian wave 
packets, in particular, the ZB phenomenon and the influ- 
ence of the initial pseudospin polarization on the charac- 
teristics of trembling motion. 

i). Following Rcf. [12] let us firstly consider the model 
problem when the lower component of initial wave func- 
tion is equal to zero, i.e. the parameters ci = 1, C2 = 
in Eq.(8a). That means that at the initial moment of 
time the electron probability is located at the sites of the 
sublattice A. It is not difficult to show that this packet is 
formed by the states with positive and negative energies. 
The relative weight of these states is equal to one. The 
wave function for i > can be found using Eqs.(9), (10): 

where the functions <j)i{f,t), (t)2{r,t) are defined by 
Eqs.(13),(14). 

In Fig.l we represent the full electron density at the 
moment t — 7 for initial wave packet, Eq.(8b) with 
width d — 2 nm and fco ~ 0.6 nm~^ . As one can 
see, at t > this packet splits in two parts moving 
along y axis with opposite velocities so that the elec- 
tron probability density is symmetrical with respect to 



FIG. 2: The average coordinate x{t) versus time (ro = d/u) 
for the wave packet with initial pseudospin polarization along 
z axis for two values of a. 

y: p(x,y,t) = p{x,—y,t) (note that at the case fco = 
the electron probability density has a cylindrical sym- 
metry at all time). For enough large time the width of 
both parts of the packet increases with time due to ef- 
fect of dispersion. One can check that in this situation 
the contributions of two components of wave functions 
ipi (f , t) and 'ip2 {r, t) in full electron density are equal. In 
other words the electron probability distributes with the 
time on the sides of sublattice A and B. Note at the 
same time p{x, y, t) ^ p{—x, y, t) and the packet center 
oscillates along x direction (zitterbewegung) . 

To analyze this motion we find the average value of 
position operator. To do it, we use the momentum rep- 
resentation. The upper {Ci{p,t)) and lower (C2(p, t)) 
components of wave function (15) in this representation 
can be easily obtained from Eqs.(ll), (12). After that 
the usual definition 



m = ^ / dpc;ip,t)zh^^^, (16) 

readily yields 

y{t) = 0, (17a) 



x(t] = ^—^ e""' [ e-'^\os(2qt]Ii(2aq)dq, (17b) 

2a J 



where Ii{z) is a modified Bessel function of the first or- 
der. In Eq.(17b) the integral term represents the zitter- 
bewegung. Note that average value x(t) depends on only 
one parameter a (in the dimensionless variables). The 
obtained functions x(t) which describes the typical tran- 
sient zitterbewegung are plotted in Fig. 2. After the os- 
cillation disappears the center of the packet is displaced 
by amount which equals to the first term Eq.(15). In the 
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case when the wave packet width is large enough and the 
inequality a — dk^ ^ 1 takes place, Eq.(17b) reduces 
to2^ 

m = 5 — (18) 

la 

As it follows from Eqs.(17), (18) for given initial po- 
larization of wave packet the ZB occurs in the direc- 
tion perpendicular to the initial momentum poy = ^^Oi 
just as for bilayer graphene^^ and for the semiconductor 
quantum well in the presence of the Rashba spin-orbit 
coupling"^' One can see from Eq.(18) that the trem- 
bling motion has a transient character as it was described 
in Refs.[12, 14] and at i > 1 x{t) I /2a. We should 
notice that Eqs.(17b), (18) coincide with corresponding 
formulas of Rcf.[14]. This is because the Hamiltonian for 
the system with Rashba-coupling 




FIG. 3: (Color on line). The electron probability density 
p{r,t) — IV'iP + l^al^ for initial wave packet, Eqs.(8a), (8b) 
with ci — 1 and C2 — 1 for a = kod = 1.2 at the time t — 7 
(in the units of d/u). 
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where a is a Rashba coupling constant, transforms into 
Hamiltonian for monolayer graphene, Eq.(l), if we make 
the replacement in Eq.(19) 



-y 



y 



u, 



(20) 



ii). Let us consider now the case when ci = C2 ~ 1; 
that is pseudospin is directed along x axis at t = 0. Then 
from Eqs.(9), (10) 



^{r, t) 



^\ (bl{f,t)+(j)2{f,t) 



(21) 



Fig. 3 illustrates the corespondent electron probability 
density at the time moment t = 7 for initial wave packet, 
Eq.(8b), for the same parameters as in Fig.l. One can 
see that the initial wave packet at t > 0, as in previous 
case, splits into two parts propagating along y in opposite 
directions so that the symmetry concerning this axis, i.e. 
p{x, y, t) = p[x, —y, t), retain during the time (as the case 
i)). The distribution of the probability density along x 
axis clearly demonstrates that its maximum is displaced 
in the positive direction that corresponds to the motion 
of the packet centre along x axis. The velocity of such 
consists of both constant as well as oscil- 



motion v,r 



dt 



latory parts. Really, a straightforward calculation yields 
the average value of position operator x 
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X / e"'? sin(2qt)—h{2aq)dq, (22) 
J a,q 





FIG. 4: The average coordinate x{t) versus time (ro — d/u) 
for the wave packet with initial pseudospin polarization along 
X axis for different values of a. 



and y{t) — like for the case i). In Fig. 4 we demonstrate 
the dependence x{t) for various values of parameter a. 
When the parameter a increases, the amplitude of ZB 
and the period of oscillations decrease. At a ^ 1 we 
have from Eq.(22) 



x(t) = ^ + sin(2at). 



(23) 



We see that the character of motion of wave packet is 
changed. Now the center of wave packet moves along x 
direction with constant velocity, which is determined by 
the first term in Eqs.(22), (23) and performs the damping 
oscillations. This constant velocity has a maximum value 
which is equal to 1/2 (in the units of u) for the narrow 
- width wave packet, or for the case fco = 0. When the 
width of packet is increased the velocity of motion of its 
centre is decreased. The frequency and amplitude of the 
zitterhewegung for a ^ 1 are the same as in the case i). 
However, the first term in Eq.(22) corresponding to the 
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FIG. 5: (Color online). The electron probability density 
p(f,t) — \4>if + IV'zl^ for initial wave packet, Eqs.(8a), (8b) 
with ci = 1 and C2 — i for a — kod — 1.2 at time t — 7 {in 
the units of d/u). 



motion of wave packet with constant velocity reduces the 
effect of ZB at least for a < 1 (Fig.4). 

It is not difficult to show that as in the other two-band 
systems the phenomenon of ZB in graphene is a result 
of an interference of states corresponding to positive and 
negative eigenenergies of Hamiltonian, Eq.(l). For wide 
enough packet a = kgd 3> 1 and at time i > 1 when 
the ZB disappears two split parts of initial wave packet 
(see Fig. 3) move along y axis with opposite velocities u 
and —u. In this situation the subpackets moving in the 
positive and negative directions consist of the states with 
positive and negative energies correspondingly. 

iii). When the initial pseudospin is along y axis the 
wave function at t > has the form 

In Fig. 5 the full electron density for the same moment 
of time and for the same parameters as in previous cases 
is shown. As one can see, the initial wave packet does 
not split into two parts at t > unlike in the cases i) 
and ii). This result is confirmed by the straightforward 
calculations. Indeed, one can show that the eigenenergy 
states corresponding to propagation in the positive direc- 
tion along y axis give the dominant contribution in total 
wave function, Eq.(24). For wide packets (a 1) almost 
all of these states belong to the positive branch of energy. 

The results of calculations of average values of x and 
y for this polarization lead to 

x{t) = 0, (25) 




0-S I 1.5 3 XJ 3 

FIG. 6: The average coordinate y{t) versus time (ro = d/u) 
for the wave packet with initial pseudospin polarization along 
y axis for different values of a. 

oo 

X f e-i\m(2qt)Ii(2aq) — . (26) 

J q 



Thus in the considering case the wave packet propagates 
along y axis and the zitterbewegung has the "longitudi- 
nal" character—. 

Note that in a numerical work^^ by Thaller the authors 
have observed similar oscillatory behavior of the expecta- 
tion value of the position operator for one - dimensional 
relativistic electron in vacuum. In Ref.[18] it was also 
shown that apart from the rapid oscillation, the wave 
packet drifts slowly even when its average momentum is 
zero. 

In Fig. 6 we represent the dependence y{t) for different 
values of parameter a. As one can see, even at zero value 
of a the oscillations are absent. In this case, as it follows 
from Eq.(26), the drift velocity is equal to 1/2 (in the 
units of u). As above, Eq.(26), takes more simple form 
at a 2> 1 



y{t) = t + ^e"* sin(2at). (27) 

Comparing Eqs.(18), (23), (27), we see that the am- 
plitude for the "longitudinal" zitterbewegung is much 
smaller than the amplitude of "transverse" zitterbewe- 
gung. This fact can bee seen as a consequence of special 
form of the initial wave function, which in the given case 
consists of (at a ^ 1) the states with positive energy 
mostly. That makes the interference between the posi- 
tive and negative components difhcult, i.e. decreases the 
ZB. Moreover, at any values of the parameter a the in- 
tegral term in Eq.(26), corresponding to the oscillating 
motion, is negligible in comparison with the first term, 
and one may neglect the effect of the "longitudinal" ZB. 

As was demonstrated above, the direction of the aver- 
age velocity depends not only on module of the compo- 
nents ipi^fyO) and ip2{r,0), but also on their phases. In 
general case for the initial Gaussian packet 
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FIG. 7: (Color online). The electron probability density 
p{f,t) = l^/'ip + |'i/'2p for initial wave packet, Eq.(8a), (8b) 
with ci = 1 and C2 = e^'^^'^ for a = kod = 1.2 at the time 
t — 7 {in the units of d/u). 
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FIG. 8: The trajectories of the center of electron wave packet 
described by Eq.(28) for two initial phases <^ = 7r/4 and 37r/4. 
The parameter a = 6. 



obtain (in the dimensionless variables) at t = 0: 

^^a:(0) = COS(^, Wy(0) = sini^. (30) 

The components of the velocity for a large enough time, 
when the trembhng motion stops, can be found from 
Eqs.(22),(26) and (29) for arbitrary parameter a 

_ 1 — exp(— a^) _ /, 1 exp(— a^)\ 

= 2a^ ^ 2^+^^ ) 

(31) 

In particular, as it follows from Eq.(31) for a ^ 1, the di- 
rection of the motion of wave packet center at large time 
coincides with the initial one, Eq(30). In other limiting 
case a ^ 1 (and for not too small (p) asymptotic direc- 
tion of the average velocity is along OF axis, i.e. along the 
average momentum of wave packet — Kk^. Thus, by 
changing the initial phase (/s, one can govern the packet 
motion and consequently the direction of dc current. To 
illustrate this, we plot in Fig. 8 the packet center trajecto- 
ries for two initial phases: = 7r/4 and 37r/4. Note that 
the packet motion with the constant velocity predicted 
above (see Eqs.(22), (26)) should lead to the existence of 
the direct current in the corresponding direction. 

To check our formalism let us consider the plane wave 
as the starting point. In this case it is easy to obtain the 
expression for the average value of electron velocity v(t). 
Really, in the Heisenberg picture the kinetic velocity is 
(in the dimensional variables) 

^^(i) = ^[f,i7] =Ma(i), (32) 

where 

do 1 ,^ -V, 2it-^ 

— = -[a,i/] = — [pxa]. (33) 

dt ih n 

In these equations p(i) — p{0). Let the initial momen- 
tum poy — hkQ. Then, using the solutions of Eqs.(32), 
(33) we find 



^{f^,0) = l^Q^y (28) 

the probability density becomes asymmetric and the av- 
erage position operator has two components 

f (t) = x{t) cos (p Ox + y{t) sin (f By, (29) 

where ip is an arbitrary phase difference between the up 
and low components of wave function and x(t), y(t) are 
determined by Eqs.(22), (26). For illustration we show 
in Fig. 7 the electron probability density obtained for the 
initial packet, Eq.(28), with (p = 7r/4. 

It is clear that the phase p determines the direction 
of the average velocity of the packet center. Using the 
expression for velocity operator v = ua and Eq.(28) we 



Vx (t) = uaz (0) sin tot + ua^ (0) cos cut, (34a) 



vy{t)^uay{0), (346) 

where lu — 2ukQ and (7^(0) = (7^ - Pauli matrixes (j = 
1,2,3). So, if in the initial state pseudospin is along z 
direction, i.e. f7z(0) — 1 (case i)) we obtain from Eq.(34a) 
that Vx (t) — u sin tot which leads to 

11 

x{t) = const cosLot. (35) 

UJ 

Returning to the original variables in Eq.(18) and setting 
d = oo we see that this expression coincides with Eq.(35). 
We get similar results also for other initial polarizations. 
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IV. CONCLUDING REMARKS 

We have studied the quantum dynamics of charge 
particles represented by Gaussian wave packets in two- 
dimensional single layer of carbon atoms (graphene). We 
investigated numerically also the spatial evolution of the 
initial wave packet and demonstrated the effect of the 
packet splitting for the pseudospin polarization perpen- 
dicular to the average momentum. The analytical expres- 
sions for the average values of position operators were ob- 
tained. These expressions clearly demonstrate that the 
evolution of wave function is accompanied by the zitter- 
bewegung and strongly depends on the initial pseudospin 
polarization. In particular, if the initial pseudospin po- 
larization coincides with initial average momentum, the 
packet center moves and exhibits the ZB along the same 
direction. In this case the second term in Eq.(26) describ- 
ing the longitudinal oscillations (the "longitudinal" ZB) 
is essentially smaller than the first one connected with the 
motion with constant velocity. As for other systems with 
two-band structur o^^i^"^'^^ , the ZB in monolayer graphene 
has a transient character. 

It was also predicted that apart from the packet center 
exhibits the trembling motion it can move with constant 
velocity (for the pseudospin polarization along x and y 
axis) . The direction of this velocity depends on not only 
the orientation of average momentum po, but also on 
the module of the components 'ijji{f,Q), tp2{r,0) and the 
differences of their phases (see Eqs.(28),(31)). 

All above calculations have been done for the K point 
of the Brillouin zone in graphene. Similar results can 



be found for initial wave packet with wave vector k in 
the valley centered in inequivalent point K'. The Dirac 
Hamiltonian around K' point can be written as^l 

Hk^=u( „ . ~P^-^Py], (36) 

This expression can be obtained from Hamiltronian 
around K point given by Eq.(l) by replacement px — > 
~Px- Thus values x{t) for the wave packet of different 
polarizations (and corresponding components of velocity) 
change sign while y{t) remain unchanged (see also^). 

In conclusion we would like to stress that the packet 
motion with the constant velocity (see Eqs.(22), (26)) 
leads to the appearance of the dc current. For the ex- 
perimental detection of this current one needs sensitive 
current meters. Experimental observation of trembling 
motion is currently a more difficult task since it is neces- 
sary to use femtosecond techniques i^ii^ If new methods of 
formation of wave packets with different pseudospin po- 
larizations will be elaborated then their trajectories and 
spatial separations can be observed probably with the 
help of devices with quantum point contacts and gates 
(see for example^^). 
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